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Control of transport barrier relaxations by resonant magnetic perturbations 
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Transport barrier relaxation oscillations, in presence of static Resonant Magnetic Perturbations 
(RMPs), are investigated using 3D global fluid turbulence simulations at the edge of a tokamak. It 
is shown that RMPs have a stabilizing effect on these relaxation oscillations and that this effect is 
due mainly to a modification of the pressure profile linked to the presence of both residual magnetic 
island chains and a stochastic layer. 
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In magnetic fusion plasmas, turbulent transport re- 
duces the energy confinement time, leading to a low con- 
finement regime (L mode). 30 years ago, the H mode, 
a high confinement regime, linked to a sheared rotation 
of the plasma, was discovered. In this regime, the turbu- 
lence is locally reduced within a narrow layer, located a 
the plasma edge, known as a transport barrier. However, 
the transport barrier is not static but relaxes quasiperi- 
odically, leading to violent bursts of heat and particles 
radially outward termed Edge Localized Modes (ELMs). 
The operationnal regime for next step fusion devices such 
as ITER will be Elmy H mode. ELMs are beneficial for 
power exhaust, but they represent however a threat for 
the plasma facing components and therefore need to be 
controlled. Over the last decade, the possibility of con- 
trolling ELMs has become more and more plausible, as 
recent experiments were carried out on DIII-D using I- 
coils, on JET using error field correction coils and on 
the TEXTOR tokamak using the dynamic ergodic di- 
vertor [3, H, 0, Q ■ These experimental studies obtained 
a qualitative control over the ELMs by imposing static 
Resonant Magnetic perturbations (RMPs) at the plasma 
edge. However, in order to get any quantitative result, 
work has to be done in the understanding of the inter- 
play between transport barrier relaxations and RMPs. 
Recently, elmy-like relaxation oscillations of transport 
barriers have been observed in global turbulence simu- 
lations of the tokamak edge [1, 0, H, 0] ■ The relaxation 
oscillations observed in these simulations have common 
characteristics with type-Ill ELMs. In this letter, we in- 
vestigate numerically the effects of static Resonant Mag- 
netic Perturbations (RMPs) on transport barrier relax- 
ation oscillations. The simulation domain is located in 
the edge region, around the g = 3 resonant surface. We 
use the geometry of the TEXTOR tokamak, and plasma 
parameters close to those used in typical experiments on 
this machine [lO, uM ■ 



The model used in this study describes electrostatic 
resistive ballooning turbulence involving the pressure p 
and the electrostatic potential 4>. The equations describ- 



ing this model are the following 
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dp 



dt 



+ {0,Vi0}= -Gp -V|(/. + z.Vi0 (1) 
{0,p} = <5eG0 -fxiiVfp +x±Vip + 5 (2) 



The first equation corresponds to the vorticity equa- 
tion, where is the vorticity of the perpendicular (to 
the magnetic field) component of the E x B flow, and 
the parallel current and viscosity effects (v) are taken 
into account. The second equation corresponds to en- 
ergy conservation, where x\\ ^'^d x± are collisional heat 
diffusivities parallel and perpendicular to the magnetic 
field, and S — S{x) is an energy source modeling a con- 
stant heat flux density from the plasma core. Following 
the standard convention, x represents the local radial co- 
ordinate, y is the local poloidal coordinate and z is the 
local toroidal coordinate, in a magnetic fusion device. 
The curvature operator G arises from the toroidal geom- 
etry of the tokamak, and Sc — 1 is basically the 
ratio of the pressure gradient length Lp to the tokamak 
major radius Rq. In the present model, time is normal- 
ized to the interchange time Tinter — c^^' ' ^"^'^ 



which 



also defines the perpendicular length scale through the 

where c, is the 



Bo 



ballooning length ^^ai = 
plasma sound speed, _Bo is the magnetic field strength, 
n\''^ is a reference ion density, mi is the ion mass, 77 1| is 
the parallel resistivity of the plasma and is a reference 
magnetic shear length. 

In order to obtain a transport barrier, the poloidally and 
toroidally averaged component 'V{x,t) = ^ Gj, of the 
E X B flow is driven by a forced poloidal sheared flow 
Vi? (a;) — ^Id tanh(x /d) By, chosen to be centered at 2: = 
(e.g. at the position of the main resonance of the RMPs), 
where fl denotes the shear rate and d is the shear-layer 
width. Here, 0(a;, t) = J^^ J^^ </> dydz is the poloidally 
and toroidally averaged electric potential. 
Note that the parallel gradient is V|| = V||o+{^_rmp , •}, 
where V||o is the component due to the unperturbed 

de- 
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magnetic field, {Vi^MP, •} = '-^Wy " ^^aS 
note the Poisson brackets. The magnetic flux due 
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to the resonant magnetic perturbation is written as: 
■>pRMpix,y,z) = j§TYl,^m{x)cos(^^^y - j^z^ where 

i^mix) = (jt) *^XP (^2;^ is the spectrum of 

the RMPs in slab geometry, mo is the central poloidal 
harmonic number, Vc denotes the radial position of the 
RMP-producing coils, is a reference value of the cur- 
rent in the RMP-producing coils, no is the toroidal har- 
monic number and rp is a typical radius where the turbu- 
lence considered in this paper (resistive ballooning) de- 
velops. In the case studied here, we use uq — 4, with 
qx=o = 90 = 3, so that the central poloidal harmonic 
number is mp — 12. 

Let the pressure be decomposed into a mean part p 
and harmonics p ^ p~ p, where p{x, t) — Jq^ Jq^ P dydz 
is the poloidally and toroidally averaged pressure. In a 
steady-state dp/dt ~ 0, the conservation of energy in 
the plasma where the conserved total energy flux is 
Qtot = J S{x)dx = const., leads to the following equa- 
tion: 

Qconv + Qdiff + QrMP = Qtot (3) 

where Qrmp — X\\{—-^-^-—'^\\P) represents the heat 
dy 

flux due to the magnetic flutter generated by the RMPs, 
and Qconv — {pvx) is the timc-averaged convective heat 
flux, where denotes fluctuations of the radial velocity, 
and Qdiff — —X-Ldp/dx is the diffusive heat flux. Here, 
(. . .) denotes an average over the poloidal direction (y), 
the toroidal direction (z), and time t. 

Figure [T] shows time series of the energy content J p dx 
and the convective flux Qconv > in presence of an imposed 
mean sheared flow, for different values of the divertor 
current Id. 

In the reference case without RMPs [Fi^T^] , by impos- 
ing a forced sheared flow on the system, a steady- 
state is reached, where so-called relaxation oscillations 
are observed, corresponding to quasi-periodic relaxations 
of the pressure gradient. These relaxations are syn- 
chronous to the heat bursts observed on the heat-flux 
time series [Fig. [T}d] and therefore also correspond to re- 
laxations of the transport barrier. 

In the case with RMPs, the energy content shows that 
the relaxation oscillations are suppressed by the RMPs, 
and this suppression is more efficient for higher values 
of the perturbation current Id [Fig. [T}:[T]3]. This sup- 
pression of relaxations is also shown as a reduction of 
the amplitude of the heat bursts and an increase of their 
frequency [Fig. [l]i[Tr]. 

Figure [2] shows the radial profiles of the pressure gradi- 
ent 1^1, the convective flux Qconv and the RMP-induced 
flux Qrmp, in presence of an imposed mean sheared flow, 
for different values of the divertor current Id ■ In the ref- 
erence case without RMPs {Id = 0) but with a mean 
shear flow (17 = 4) , the convective flux Qconv [Fig. [2)3] is 
reduced around the position a; = compared with a case 
with no flow [Fig. [2)3, dash-dotted line]. A high pres- 
sure gradient |^| around this position, corresponding to 




FIG. 1: Effects of RMPs on the dynamics of barrier relax- 
ations: time traces of the thermal energy (a,c,e) and the ra- 
dial heat flux (b,d,f), in presence of an imposed mean shear 
flow with shear rate fl = 4, for (a, b) no RMP perturbation, 
(c, d) In = 0.5 kA and (e, f) Id ^ 1 kA. 

an ETB is created at the position of maximal flow-shear 
X — Q [Fig. [2^] . The appearance of this strong pressure 
gradient is linked to a reduction of the convective heat 
flux Qconv by the mean sheared flow as seen from the 
conservation of energy ([3]) with Qrmp = 0: 

Q^conv + Q^diff = Qtot (4) 

where the superscript indicates the reference case 
without RMPs {Id = 0). 

In the case when there is a combination of a mean shear 
flow {n = 4) and RMPs {Id ^ 0), the ETB generated by 
the mean shear flow is eroded in the vicinity of the posi- 
tion re ~ 0, compared to the case without RMPs [Fig. [2^]. 
This erosion of the ETB only appears when both a mean 
shear flow and RMPs are present, and therefore can be 
explained by a synergetic effect, p 

We propose the following model based on the balance 
of heat fluxes, to explain the behaviour of the convec- 
tive flux and pressure gradient in the presence of RMPs 
and a shearflow-induced transport barrier. Taking into 
account RMPs, it is convenient to decompose the pres- 
sure harmonics into an equilibrium and a turbulent part: 
p = p'^'^{x, y, z) +p*'^^^{x, y, z, t), and similarly for the ra- 
dial velocity harmonics Vx- The energy balance ([3]) can 
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\J Id / depending on the region considered, there are 
three different hmits: 

i) Close to the resonant surface r — tq, e.g. ^ 
^ ^ d, the turbulent convective flux is reduced by the 
sheared flow and additionnally reduced by RMPs due 
to stochasticity of magnetic field lines [Fig. [4] so that 
Qmm/Qtot ^ 1 [Fig. [3K]. Moreover, an equihbrium 
convective flux Qlonv indirectly linked to RMPs appears 
[Fig. [31d]. This equilibrium convective flux is a conse- 
quence of the presence of an RMP-induced residual mag- 
netic island chain at the radius r = tq (in TEXTOR, 
ro = 0.45 [m]) corresponding to the position x = 
[Fig. [4]. Moreover, the RMP induced flux is low in the 
vicinity of x = 0: QBMp/Qtot < 1 [Fig. [2t]. Therefore, 
the energy balance ^ simplifies to: 
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FIG. 2: Effects of resonant magnetic perturbations on: (a) the 
pressure gradient and (b,c) radial heat fluxes, for different 
values of the perturbation current Id- Qconv denotes the 
convective heat flux and Qrmp is the conductive heat flux 
induced by the resonant magnetic perturbations. The total 
heat flux is Qtot = 10. 
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FIG. 3: Components of the radial convective flux: (a) tur- 
bulent convective flux Q^^^v, and (b) equihbrium convective 
flux Qtonv, e.g. RMP-linked, for different values of the per- 
turbation current Id- The total heat flux is Qtot = 10. 



so that, in the region a; 0, the appearance of the equi- 
librium convective flux Q^lnv must be balanced by a de- 
crease of the pressure gradient, seen on Figure [2^, similar 
to a flattening of the pressure proflle on the island chain, 
that occurs in the study of plasma macroinstabilities such 
as tearing modes [l^ . 

ii) For r < rg, far from the resonant surface but inside 
the barrier region, e.g. —d <^ x <^ — ^, the RMP-linked 
equilibrium convective flux is small Q'^^nvlQtot <SC 1, 
since there is no residual island chain in this region, 
so it does not play any role. The turbulent convective 
flux Q*"^^ is reduced by the mean-shear flow. More- 
over, the RMP-induced flux Qrmp is small in this re- 
gion Qrmp /Qtot ^ 1. Therefore, the energy balance ^ 
reduces to: 



-\turh 
^ conv 



idiff 



(7) 



Additionnally, in this region, the turbulent convective 
flux does not depend on the perturbation current Id 
[Fig. [3^] because the magnetic field lines are not stochas- 
tic [Fig. m, so that for <C a; < the turbulent 
convective flux is only weakly perturbed: Qlo^^ ~ Qconv 
Therefore, equations (g]) and ([7]) imply: Qdiff ~ Qdiff 
Thus the pressure gradient profile is only weakly modi- 
fied by the RMPs [Fig. Hi]. 

iii) For r > rg, far from the resonant surface but inside 
the barrier region, e.g. ^ ^ a; ^ d, the equilibrium con- 
vective flux is small {Qtlnv/Qtot ^ 1), since there are 
no residual island chains. The energy balance ([3|) thus 
reduces to: 



then be written: 



^ conv 



-\turh 
-> conv 



Uiff ■ 



Qrmp = Qt 



(5) 



where QH^^ 



ip'^'^v^^) is the equilibrium convective flux 
) is the time averaged convective 



~turb~turb 



and QlZl = {P- 
turbulent flux. 

Because of the presence of localized residual island 
chains (as seen on Figure S]) of normalized width w ~ 



QlTnv + Qdiff + Qrmp ^ Qt 



(8) 



An increase of the perturbation current Id causes an 
increase of the RMP-induced heat fiux Qrmp, but it 
also induces a decrease of the turbulent convective heat 
flux linked to stochasticity of magnetic held lines 

[Fig. g]. In this region, there is therefore a competition 
between the Qrmp and beat fluxes, which may 

explain the fact that, in the region ^ <C cc <SC d, the pres- 
sure gradient increases for small perturbation currents 
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FIG. 4: Effects of RMPs on the topology of magnetic field 
lines: Poincare map for a perturbation current Id ~ 0.5 kA. 
There is clear evidence of three residual magnetic island 
chains, and of a stochastic layer in between the central is- 
land chain and the one on the right. 
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FIG. 5: Effects of width d of the shear-layer on the dynamics 
of barrier relaxations: time traces of the radial turbulent flux, 
in presence of an imposed mean sheared flow with shear rate 
n = 2, for a) d = 20% and 6) d = 10%. 



and decreases for higher perturbation currents [Fig. [2^]. 

From the analysis presented here, the stabilization of 
relaxation oscillations appears to be linked to an erosion 
of the transport barrier, i.e. a reduction of its width. 

To confort this hypothesis, we performed simulations 
without RMPs, but with a narrower transport barrier 
(produced by a sheared flow with a thinner shear-layer d) . 
Figure [5] shows time series of the convective flux, in the 
case without RMPs, for two different values of the shear- 
layer width d. As seen from the comparison of Fig. [1] and 
Fig. m there is a striking similarity between the effects of 
a decrease of the shear-layer width d and the effects of 
RMPs on the dynamics of Edge Transport Barrier (ETB) 
relaxations. Our simulations therefore suggests that the 
main effect of RMPs in presence of an ETB is to modify 
the geometrical properties of the ETB, e.g its width and 
position, yielding a reduction of the amplitude and fre- 
quency of the elmy-like relaxations, and therefore leading 
to grassy elmy-like relaxations. 

To summarize, we investigated the effects of resonant 
magnetic perturbations (RMPs) on transport barrier re- 
laxations. In conclusion, it is shown that RMPs have 
a stabilizing effect on these relaxations. This effect is 
linked to a modification of the pressure gradient equilib- 
rium profile due mainly to a modification of the magnetic 
topology, e.g. the formation of magnetic island chains, 
and to RMP-induced stochastic transport. An erosion of 
the pressure gradient profile is observed at the surface of 
principal resonance (also chosen to be the central posi- 
tion of the sheared flow). This erosion is shown to be 
linked to the presence of residual magnetic island chains 
inducing a stationnary convective transport of heat (and 
particles) in the radial direction. Far from the principal 
resonance surface but inside the shear-layer, the pressure 
gradient modifications are only linked to the presence (or 
not) of stochastic resonance overlap. 
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